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for  this  purpose.  The  theory  requires  the  calculation  of  the  coefficient 
of  a certain  quadratic  form  in  ei^t  independent  variables.  Explicit 
formulas  are  given  for  the  most  complicated  of  these  coefficients.  They 
are  indeed  so  complex  that  the  author  has  not  succeeded  in  drawing 
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ble for  the  numerical  analysis  of  particular  cases. 
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1.  Introduction.  In  two  previous  papers,^  hereafter  denoted  by 
and  MP2,  the  author  has  developed  a "metric  theory"  of  ordinary 

differential  eq.uations  with  the  view  of  studying  the  dependence  of 
the  solutions  on  the  initial  conditions.  In  Section  4 of  M?2  some 

indications  were  given  as  to  how  the  theory  could  be  applied  to  the 
eiiuations  of  exterior  ballistics.  The  purpose  of  this  paper  is  to  go 
into  more  details  with  regard  to  this  application.  The  ultimate  goal 
is  to  obtain  results  suitable  for  numerical  analysis  as  well  as  for  a 
clear  iUalitative  picture  of  what  can  be  expected  in  problems  of  this 
sort. 

2.  The  eciuations  of  motion  of  a spinning  symmetric  shell.  We  use  a 
moving  cartesian  reference  ftrame  F with  origin  0 at  the  center  of  gravity 
of  the  shell  and  with  one  of  its  axes  (namely  the  XQ-axis)  coinciding 
with  the  shell  axis  and  directed  from  0 to  the  shell  vertex.  Ihe  other 
two  axes  (namely  the  xi-axis  and  X2-axis)  are  chosen  in  such  a way  that 

0-XjX^Q  form  a right  handed  system.  Tt  is  not  necessary  to,  assume  that 

the  shell  is  rigidly  attached  to  this  reference  frame  F.  In  fact  we 
slaqpllfy  owe  equations  slightly  if  we  assume,  as  we  do  frem  now  on,  that 
the  component  of  the  angular  velocity  of  F in  the  direction  Oxq  is  zero, 2 
while  that  of  the  shell  is  cOq,  in  general  / 0.  On  the  other  hand  the 

components  of  angular  velocity  of  both  F and  the  shell  in  the  direction 
of  the  other  two  axes  are  the  same  and  wiLl  be  denoted  by  u)^  and 

respectively.  Hence,  if  is  the  vector  angular  velocity  of  F 

(relative  to  some  fixed  frame  Fq)  and  if  is  the  vector  angular 

velocity  of  the  shell,  we  have 

A 

” Oi 

vhere  i,  J,  k are  unit  vectors  directed  along  the  Xr>,  x,,  x„  axes 
respectively.  ^ 1 2 

Similarly  we  denote  the  vector  velocity  V of  the  point  0(relative 
to  the  fixed  frame  Fq)  la  the  form 

V = ’JqI  + Ugk, 

while  the  vector  G representing  the  acceleration  of  gravity  and  having 
constant  components  relative  to  the  fixed  frame  Fq  appears  in  the  form 

G “ Kq!  + 4-  ggk. 
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Denoting  the  mass  of  the  shell  by  m,  its  axial  moment  of  inertia  by  a, 
and  its  transverse  moment  of  inertia  by  b,  we  can  write  the  moment  of 
momentum  in  the  form, 

Jl/J^  - a + b t b o^gk 

and  the  momentum  appears  in  the  form  mV.  We  denote  the  vector  aerodynamic 
moment  by 

M = aj^i  + bMj_J  ■+■  bMgk, 

the  factors  a and  b being  introduced  for  later  convenience;  and  we  denote 
vetior  aerodynamic  force  by 

m ^ = m -f  m^^j  + ^^2^’ 

We  now  are  in  a position  to  write  down  three  vector  differential 
equations.  The  first  of  these  is 

(2.1)  mV  f + m V = mG  + m ^ , 

which  e.iuates  the  rate  of  change  of  momentum  (relative  to  the  fixed  frame 
FQ)to  the  resultant  of  the  applied  forces.  The  second  of  our  e^iuatlons  is 

M ^ , 

which  eciuates  the  rate  of  change  of  moment  of  momentum  to  the  applied 
tor^iue.  Finally,  the  third  e.^uation 

(2.3)  6 i “ +i2p>f  0 • 0 

merely  expresses  the  fact  that  G 15  a constant  vector;  that  is,  it  has 
constant  components  relative  to  the  fixed  frame  Fq. 

Now  ^ and  M are  assumed  to  be  functions  of  J2^  and  V only,  !Die 

precise  n.iture  of  this  dependence  will  be  specified  in  detail  later  on. 

We  only  wish  to  remark  now  that  the  above  three  vector  e>iuations  become 
self  contained.  They  are  equivalent  to  nine  scalar  differential  equations 
for  the  determination  of  nine  unknown  functions,  namely  0^2^/ a;g> 

Ug,  gQ,  gj^,  gg.  Actually  the  order  of  the  system  can  be  reduced 

from  nine  to  eight,  because  of  the  obvious  first  integral  G • G - conet. 
When  once  these  unknown  functions  have  been  determined  by  integration 
of  (2.1),  (2.2)  and  (2.3),  the  trajectory  may  be  calculated  as  follows: 

We  introduce  a moving  frame  Fq'  with  the  same  origin  as  F (namely 

the  center  of  gravity  of  the  shell)  but  with  its  axes  invariably  parallel 
to  the  corresponding  axes  of  the  fixed  frame  Fq.  Eulerian  angles  Q,  (p  , y/ 

are  now  Introduced  to  define  the  orientation  of  F relative  to  Fq'  and 
hence  also  relative  to  Fq.  We  hereby  assume  that  the  axes  OJ,  0^,  oy 
of  Fq'  form  a right  handed  system,  that  9 is  the  angle  from  O^to  OXq, 
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that  is  an  angle  ftron  Oj  to  the  intersection  ON  of  the  -plane 
with  the  Xj^Xg-plane  (the  so-called  line  of  nodes),  and  that  is  an 
angle  from  this  scoae  line  of  nodes  ON  to  the  Qx^  axis.  Cf.  the  figure. 


We  further  assixme  that  F_ 
and  hence  F^'  is  chosen  in  su& 

a way  that  one  of  its  axes,  say 
the  i -axis,  points  vertically 

upward.  !Ihen,  since  gravity 
always  acts  vertically  downward 
we  see  easily  that, 

g^  = - |G|  cos  0 

® siny' 
gg  = - /G|  sin  0 cosy 


/ 2 2 

from  which  we  find  that  tan  y - gj^/gg  and  tan  0 ■ -f  j/g^  4-  /g^. 

When  yr  and  0 are  thus  found  as  functions  of  t (since  gQ,  g^,  g^  are 


known  functions  of  t after  the  integration  of  (2.1),  (2.2),  and  (2.3)), 
we  ceui  find  <p  by  the  single  laadrature  Involved  in  integrating  the 
equation 


^ ^ cos  e « 0. 

3t  dt 


This  equation  expresses  the  fact  that  the  frame  F has  zero  nTigniay 
velocity  about  Oxn.  Thus 

t 

(2»^)  <p  - ~ ^ Y ® 

■4 

With  the  Eulerian  angle  all  determined,  we  can  now  find  the  nine  direction 
cosines  between  the  frames  F and  Fq'  and  hence  between  F and  Fq.  They 

are  exhibited  in  the  following  table. 


X^^cosfcoaf-  slnfslnjrcog  eU2-»ln^co*)r+  cos f Bln)/ cos  6 ^sslnjf'sln  0 


»±n.yr±  cos ^ cosy  cos 


XQlV^j^ssin  e Bin  <f> 


Vg^-sln  e coB<p 


yycos  0 


vfaCM  the  axes  of  the  fixed  ftrosae  Tq  are  denoted  by  X,  T,  Z.  This  table 
my  easily  be  obtained  fron  the  lav  of  cosines  in  spherical  triconeotetry 
or  by  the  aatrix  Kiltlpllcation  corresponding  to  the  fact  that  F nay  be 

obtained  frcm  F^'  by  rotating  the  latter  through  ^ about  the  X axis, 

then  throu^  0 about  the  carried  ^ axis,  and  finally  throu^  y about 

the  carried  -axis. 

ceavonents  ^ and  ^ of  the  shell  velocity  relative  to  F^ 
dt  dt  dt 

are,  of  course,  related  to  the  coaponents  u^,  xi2,  Uq  relative  to  F by 
aeans  of  the  foUovlng  equations: 


(2-5) 


i = ''a“i  * ■*  *'2«0 

i • * *'3"0- 


Since  the  direction  cosines  denoted  by  the  ^'s,  f*.'B,  and  )f 's  are  known 
functions  of  t as  well  as  the  u's  (when  once  equations  (2.1),  (2.2), 
(2.3),  and  (2.4)  have  been  integrated),  the  actual  trajectory  nay  be 
obtained  by  the  three  obvious  quadratures  poresented  in  (2.$). 

Thus  the  key  to  the  entire  situation  is  the  integration  of  the  self 
contained  system  of  ninth  order  represented  by  the  three  vector  equations 
(2.1),  (2.2),  and  (2.3).  It  is  this  system  to  which  ve  shall  apply  the 
metric  theo^  of  our  previous  papers,  b^rltten  out  in  full,  these  nine 
equations  take  the  following  form 
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> 

(2.6) 

du 

— -“'A 

dt 

(2.-7) 

^^1  4 CO  U a g J.  ^ 

ar*  20  ^i^5^i 

(2.8) 

^“2  -co  u = g + 5 
3t“  1 0 ®2  + iP  2 

(2.9) 

dt  ^ 

(2.10) 

+ h 6U  io  = M (^ere  h - a/b) 

dT  021 

(2.11) 

d^ 

— £ - hco  o>  s U 
dt  0 12 

(2.12) 

— - -fCOg  -<JOg  rO 

dt  1 2 2^1 

• 

(2.13) 

^1 

.fOOg  = 0 

dt  2 0 

(2.11^) 

— £ - to  g =0. 
dt  10 

If  we  introduce  complex  quantities  defined  as  follows,  the  number 
of  our  enuations  can  be  reduced  from  nine  to  six.  Since  three  of  these 
six  equations  are  complex,  this  simplication  is  formal  rather  than  actual. 
Still  it  is  veil  worth  carrying  out.  We  set 

u s u + iu^ 

1 2 

(2.15)  ^ (not  to  be  confused  with  vector^) 

M « + IMg 

g ■ 8X  + igg 
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Then  multiplying  (2.8)  by  i and  adding  to  (2.7)  we  obtain,  with  the  help 
of  (2.15),  the  following  csmplex  equation 

(2.16)  du  z iu  60  -f  g + 5 

dt  0 ^ 

to  be  used  as  the  equivalent  of  the  two  equations  (2.7)  and  (2.8).  In 
a similar  manner  (2.10)  and  (2.11)  are  replaced  by 

(2.17)  ^rihcodO+M, 

dt  ^ 

while  (2.13)  and  (2.14)  are  replaced  by 

(2.18)  ^ = ig  60 

dt  ^ 

3.  Expressions  for  the  aerodynamic  force  and  torque.  The  use  of  con5>lex 

quantities  also  simplifies  the  discussion  of  the  effect  of  axial  symmetry^ 
on  ^and  M«  Thus  M - uq,  u).  A rotation  of  the  frame  F through 

an  angle  oc.  about  Qxq  replaces  M by  >^by  , and  u by  ue^*^.  Hence, 

by  axial  symmetry,  we  must  have 

(3.1)  Uq,  u)e^*^=  M(o>q,  u^,  uei®^) 

as  a functional  equation  to  be  satisfied  by  M for  all  real  values  of  ot . 

If  we  consider  only  terms  that  are  linear  in  the'  relatively  small  quantities 
Ui,  andcOg,  it  is  readily  shown  from  (3*l)  that  the  most  general 

foiTD'  for  M would  be 

(3.2)  M = Au)-h  Bu 

where  A and  B are  complex  valued  functions  of  the  real  quantities  u^  and 
fjQ  . A similar  argument  shows  that 

(3»3)  ^ r C6iJ  ^ Du, 

where  C and  D are  also  conq)lex  valued  functions  of  u^  and  oO^,  and  that 

(3.^)  Mo  = Mo  ( ^<^0,  Uq) 

(3.5)  (a>0>  “o> 

where  andj^O  real  val>\ed  functions  of  xiq  and6*>o» 
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If  ve  as  SUM  that  A caa  be  develeped  In  powers  of  and4x>Q  and 
if  we  then  neglect  all  terms  of  hl^r  degree  than  the  first,  we  have 

A s (a^^  ■+  a-n  o)q  ■+  ■#'  ^(^2  ^pi  ^0  ^22^^^* 

ISivis  there  are  six  real  "aerodynamic  constants"  associated  with  A, 
namely  a^^,  ag>  a^,  a^^,  &21,  &22‘  Similarly  we  write 

B = (bj^  -f  ^u.^0  -f  ^^^2  + ^21  "^0  ^22'"o) 

C « {c^-h  <^n^0  + Ci2'^)  + i(c2  + + ®22'^^ 

D a (d^  + dii^O  **•  ^22^^) 

Mq  « Bq  ^ ”11 “12'*^0  '\)  “22'"0^ 

3o  ^ ^0  + ^1^0  ^2'^  ^11^0^  ^12  ^ ^22^^)^ 

Thus  we  consider  at  the  outset  a total  of  3^  aerodynamic  constants. 
However,  it  is  physically  clear  that  1^  ■ 0 if  to^  - 0.  Hence 

"b  " *2  “ *22  ■ °‘  ^*°^0  = 0 ^ «0  ‘ ^0  ‘ ^1  " ^11  ’ 

Also  is  an  odd  function  ofc<>Q  while  is  an  even  function  ofu)^. 

Bence  * 0 and  f^  « 0.  The  nuBiher  of  constants  is  thus  reduced  at 

once  to  28.  furthermere,  it  seems  to  be  generally  agreed  that  of  these 
28  constants  aj^,  a2,  bj^,  bg#  Cg,  dj^,  d^,  m^^,  and  fg  may  be  neglected. 

We  are  thus  left  with  I8  constants  and  the  above  equations  reduce  to 
A = ^ 12'*0^  + ^(^21^' 0 ®22^) 

B « '"12'*0^  •*’  ^22^*0^ 

C • (c^«>Q  ■#•  c^Uq)  f KCg^^^O  + ‘^22"'0^ 

D = (*1x^0  ^22'^0^ 

**0  ‘ “12^0^^0 

4 ' ^22"0^ 
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The  elimination  of  eight  more  constants  can  be  accooplisbed  as 
foUovs:  If  « VI2  B 0 equations  (3«2)  and  (3*3)  become 

M B At  -f  §•  Ci  <^2  that,  v^on  separating  real  and 

pure  imaginary  parte,  we  have  (with  the  help  of  the  above  ejqpressions 
for  A,  B,  C,  D)  that  foUeving  relations: 

*22“o>“’a  + ('’u‘^0  + Vo\ 

= (*u“>0  + •l2"o)“’2  + (’>2l‘^0  + *’22"o\ 

h ■ + 'l2"o>“'2  + <*2l“^0  '*22"0>^- 


But , with  « Ug  s 0,  the  vector  velocity  of  every  point  on  the  shell 
axis  must  be  orthegonal  to  the  X2>axis.  Hence  the  ceogponMit  in  the 
direction  of  the  xj^-axis  of  any  element  of  aerodynamic  force  acting  on 

the  shell  would  not  chanse  sign  if  the  direction  of  spin  were  reversed, 
but  the  cesqpoiient  in  the  direction  of  the  Xg-axis  woeld  change  its  sign 


under  such  circumstances.  In  reaching  this  conelusion  we  use  consider- 
atiens  of  syawtry  in  the  form  seawtlmes  called  "the  jvlnciple  of 
sufficient  reas<»"  (cf . a.  0.  Blrkhoff,  Collected  Mathematical  lepers, 


volume  3#  PP«  778-804),  We  thus  conclude  that,  with«>^  " '*2  * ®'-^l 


and  3Ktst  be  even  fUmetlons  of  o>q,  while  j^g  and  must  be  odd 


functions  of  We  thus  read  off  fresi  the  above  equations  that 


®22  " '*12  * ®^U  * ^21  " ®21  * ^11  * ®12  “ ^ final 

formulas  fOL'  A,  B,  C,  B,  1^,  and^Q  are  as  follows: 


A = * i»21% 

B s bjj^^Q  ibggUQ 


C = i«22'*0 

D * d^gUf)  *f  i^21^0 

•S,  • 
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Qhe  ten  remaining  constants  are  usually  expressed  in  terms  of 
dimensionless  quantities  i^ich  are  generally  denoted  by  the  capital 
letter  K with  various  subscripts.  In  accordance  witoi  one  well  known 


usage  we 

write 

^12 

= - 

^21  = 

'‘n 

= - 

b = 
22 

- p^V'^ 

°22  ■ 

3 

pd  Kgffi  1 

^12 

= - 

^21  ’ 

pd3K^m-l 

\2 

. - 

^22  " 

- m”^ 

r DA 

Here  (>-  density  of  air  (=  .001188  gm/caJ  for  a typical  value) 

d = diameter  of  projectile  (2  centimeters  in  a typical  case) 
m = mass  of  projectile  (175  grams) 
b = transverse  moment  of  inertia  (8OO  gm  - cm^) 

a = axial  mement  of  inertia  (76  gm  - cm^) 

According  to  some  experimental  work  of  Turetsky,  the  dimensionless 
aerodynamic  coefficients  may  be  exi^ected  to  have  values  of  about  the 
following  magnitude:  ^ ^ * -0.1,  ■ 1,  * - 10. 

Kj  ■ X,  . 0.2.  = 0.005,  = 0.1.  I have  no  informtion  about 

and  Kjjp.  Ihey  are  probably  too  small  to  be  significant. 

IJeedless  to  say,  the  aerodynamic  force  system  described  here  does 
not  satisfy  the  Nielsen-Synge  requirement  of  invariance  with  respect  to 
shift  of  mass  center. 

k.  D^esslon  on  the  approximate  integration  of  linear  differential 
equations . ftie  approximate  solution  of  the  differential  system 
introduced  above  hinges  on  the  solution  of  a system  of  linear  equations 
of  the  following  form: 


f 


>■ 


I 


h..- 


(4.1) 


~ * «(t)n'  + ^(t)u  + H(t) 

g . y(t)w'  + <^(t)v  + L(t) 
o.t 


•• 


•••  t • • «•  •• 
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vhere  Uj  fit  >(,  B>  azid  L are  conplex  valued  functions  of  the 

real  varlahle  t.  The  solution  of  euch  a system  of  non-hoaog«»eous 
linear  dlffexmntlal  equations  is  veil  known  to  depend  on  the  solution 
of  the  cerresponding  hooogeneous  system.  In  which  H(t)  and  L(t)  axe 
both  zero.  But  It  Is  possible  to  write  down  iaeedlately  an  approximate 
solution  of  the  hoieogeneous  system,  at  least  If  h^  « i^(<x  - <S’)2  + kfij 

and  are  not  zero  and  the  derivatives  of  «,  fi,  /,  S ere  small  enough 
coaqpaxed  to  h and  . In  fact  the  functions 

W = pih’t  eaq?  -h  hj  dt 

(^.2)  U = |h  - i(<x  - J )J  expj' 

may  be  verified  to  satisfy  exactly  the  equations 

f = («*if 

('‘•3)  §■  (r  -i  ±:l)w*(S^ii.iJ)n, 

E ^ E r 

which,  in  view  of  th^  smallness  of  h,  ^ and  of.--  S , are  only  slight 

ncdlficatlons  of  (4.1'  when  H(t)  « L(t)  * 0.  Moreover  fr'in  the  definition 
of  h as  a square  root’,  of  a non-vanishing  q';5jmtlty,  it  Is  clear  that  (4.2) 
yields  two  solutions,  one  for  each  de'tezmlnatien  of  the  sign  of  'the  square 
root.  fEe  verification  that  (4.2)  satisfies  (4.3)  is  laborious  but 
elementa-y  and  will  bo  left  to  the  reader. 

To  ob'tain  a solution  of  the  non-hoaogeneo;ts  linear  system  (4.1)  we 
need  to  know  two  particular  solutions  of  the  homogeneous  s^tem  depending 
on  a real  parameter  s,  namely  w^(t,s),  u^(t,s)  and  »*'2(t,8),  Ug(t,s), 

such  that  vvj^(s,s)  ■ 1,  u^(3,s)  s 0,  . 0,  Ug(B,s)  , 1.  Then,  as 

is  well  known  and  as  is  easily  verified,  it  turns  out  that  a pBirticular 
solution  of  (4.1)  is  given  by 

r 

^ " J L 4-  L(s)  ds 

(4.4) 

u - 


*^*v 


•i 


K 
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But  frm  the  apjroxlmte  solution  (4.2)  of  t]ie  honogeaeous  equations 
it  is  easy  apprOKifflately  to  set  these  functions  w'^Ct^s)  and  u^(t,s). 

m fad^  using  the  abhreriations  ^ 

P(t)  * 

Qj^(t)  « (h  - i(«-5))  /iTV^ 


Qg(t)  - (.  - h - - <T  )) 

VB  vrlte  ^ t T 

‘^(t,8)  = iP(t)  exp  ^ - ^(e)  «tp  ^ h)4*y 


(4.5.1) 


U3^(t,s)  - i j^<i^(t)q^(a)  exp^  - h)dr  - (i^(t)i^i»)exp  ^ (i(A^  -f>h)dxj 

^ t * 

v»'2(**e)  = if(t)P(s)jexp  ^ ) -h  h)dX  “ ^ 


(4.5.2) 


U2(*>®)  = iP(e)|oj^(t)eaq?  ^ (t(-< +J)  t ^)dt  - «^(t)  'itp  ^ (i(«->^^  ) - h)^-^ 

a » 

Both  vV^(t,s),  u^(t,s)  and  »v(t,s),  Ug(t,s),  as  given  by  (4.5),  satisfy 

(4.3)  exactly,  lAlle  (4.4),  with  these  Yal\»^  f«r  the  axid 
Inserted,  will  satisfy  exactly  the  foUoving  sys^ai; 


(*‘-*- - i b)iw 
ox  r h 


■^H(t) 


(4.6) 


(Y^ii  o^-r  _i  ci-M  ^ - U)u+ at), 

dt 

which,  of  course,  is  a close  appraxinetien  of  (4.1).  thus  the  functloos 
given  by  (4.4)  will  approxinate  closely  the  solution  of  (4.1)  vtdch 
vanishes  when  t ■ 0,  and  cstlnates  of  this  apjroocinatlec  can  be  carried 
out  by  classical  methods. 


■H Jt, 
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5*  The  ^tegratlon  of  the  equatl<^  of  aotlon*  We  T^in  by  writing  down 
in  a suitable  form  the  ey,uations  developed  in  Section  2 and  3.  Namely 
we  have  enuations  (2.6),  (2.9),  (2.16),  (2.17),  (2,l6),  while  (2.12)  is 

eliminated  by  setting  g * - ( )g1^  - g,?  - g ^)2.  We  also  substiti*e 

for  ^ , Mq,  and  N their  values  in  terns  of  the  aerodynamic  constants 

discussed  in  Section  3.  The  resiat  is  the  following  systea  of  five 
e.. nations: 

(5.1)  du^/dt  = - u)^u^  - (IGI  ^ - |gi 

(5.2)  dco^dt  * 

(5*3)  dg/dt  - - ltv(|G|^  - |g/^)^ 

(5.4)  dw/dt  . a^2'"o  ^^^21  ^ 

(5.5)  du/dt  = (cii^^Q  + i(c22  i.  Duq)!*'^.  (d^u^  + 8 

We  have  set  = m/  iq  order  co  conform  to  the  notation  of  the  previous 
section.  Since  the  last  three  eviviations  are  co^lex  equations,  they  are 
equivalent  to  six  real  equations.  Hence  the  order  of  our  system  (considered 
as  a real  system)  is  eight,  even  thoi;igh  for  some  purposes  it  may  be  txmated 
as  of  order  five.  Bie  quantities  l0|  , m^g,  a^,  ^2^,  h,  b^,  b^^, 

°11^  °22^  ^1^»  ^1  course,  all  constants,  while  Uq,  cOq  g ^ g^  + Ig  , 

v/z  w 2 i w)g,  ji  « u^  + iUg  are  the  unknown  functions  of  t. 


^ equations  can  be  trivially  satisfied  by  taking  g,  w/  . and 
u idMtica^  zero  and  then  by  integrating  the  first  two  equations  by 
s^le  quadratures.  »jt  this  would  correspond  to  shots  fired  either 
straight  up  or  straight  down  and  is  of  no  practical  significance. 

method  of  obtaining  exact  elementary 
solutions  of  these  equations,  although  the  classical  eSstence  theoreas^^ 
not  only  ^sert  that  solutions  corresponding  to  arbitrary  initial  values 
of  the  ^owns  do  exist  but  give  lJ5>licitly  a means  fS^the  numerical 
^cuUtion  of  these  solutions  to  any  desired  degree  of  accuracy.  These 
methods  ustol^  start  with  a.  crude  approximation  (which  we  caU^the 
approx^tion)  and  then  proceed  by  successive  approximations.  The  nti 
approximation,  namely,  is  obtained  freo  the  (n^ll^  w 

(-1)^  apprOKl«tS“lf?hi“r]it  - 

i'  then  integrating  under  the  appropriate  initial  conditions  it 

is  known  (^uider  suitable  precautions)  t£t^  approxi^tJ^STci^ige 
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unifonOy  to  the  exact  solution  and,  in  fact,  they  converge  quite 
rapidly  if  the  ti»e  interval  is  not  unreasonably  large,  no  nutter  how 
crude  the  zero^  approximation  may  be.  We  however,  focus 

attention  on  how  to  get  a really  refined  zero^  approximation,  which 
is  sufficiently  accurate  so  that  the  first  approtinution  may  be 
considered  (at  least  for  some  purposes)  as  the  actual  solution. 

We  hereby  limit  ourselves  to  initial  conditions  in  which  the  angle 
between  the  shell  axis  and  the  tangent  to  the  trajectory  is  small.  Whether 
this  so-called  "yaw-ai^e”  remains  small  depends  on  the  differential 
equations.  For  constats  correspondl^  to  a pr^rly  constructed  shell 
it  should  renuin  small,  and  our  zero^  approximation  will  in  fact  be 
obtained  by  taking  it  actually  eqiul  to  zero.  It  is  for  such  properly 
constructed  shells  that  our  system  of  successive  approximations  is  to 
be  expected  to  converge  very  rapidly.  For  lxg«ropcrly  consti  acted  shells 
the  approximations  theoretically  would  still  converge  but  in  m much  slower 
and  more  unstable  manner. 

We  denote  by  0 the  angle  between  the  tangent  to  the  trajectory  and 
the  vertical  direction  (i.e.  the  ^-axis).  If  the  yaw  angle  is  small 

then  6 & is  small  (although  the  converse  is  not  necessarily  true). 

m all  cases  Q - & is  to  be  taken  as  saull  for  t = 0,  and  in  the  zero^ 

approximation  0 for  all  t. 


A more  complete  statement  of  the  initial  conditions  for  both  the 
zerotb  approximation  and  the  actual  solution  is  indicated  in  the  following 
table: 


Variable  For  the  zeroth  approximation  For  the  actual  solution 


“0  1 

> 0 (large) 

>0  (large) 

u 

= 0 

coagplex  number  with  small  modiU.us 

COq 

> 0 (large) 

^0  (large) 

w 

• 0 

complex  number  with  small  modulus 

e 

anjlo  between  0 andTr* 

angle  between  0 and  v* 

= - |g/  cos  B 

approximately  the  same 

S2 

= - /G/  slnfl* 

approximately  the  Scime 

r 0 

small  real  nuaiber 

(p 

= 0 

small 

Y 

small 

•»thold3  also  throughout  the  motion 
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>|VtV<«g  <p , and  y equal  to  zero  initially  means  that  the 
shell  axis  is  initially  in  the  TZ  (n  y^)  plane.  If,  in  addition, 
u * 0,  the  trajectory  is  also  in  this  plane,  at  least,  initially. 

To  get  the  zeroth  approximation  ve  assume  that  the  center  of  gravity 
of  the  shell  foUovs  the  trajectory  of  a particle  in  exterior  halllstics, 
in  vhlch  the  drag  is  proportional  to  the  square  of  the  velocity,  Eind  that 
the  yav  is  zero.  This  ameunts  to  ^placing  equations  (^.l)  and  (5*3)  ^y 


(5.6) 

du 

0 

CM 

1 

■ 

- g + f u ^ 

dt 

2 22  0 

and 

(5.7) 

dt 

^0 

1 

|iS 

respectively.  Here  the  drag  coefficient  is,  of  course,  > f^g.  These 

equations  are  obtained  from  Revton's  second  lav  of  motion  by  talcing 
conpoaents  In  the  tangential  and  normal  directions;  and  (5  *7}  ^Amn 
expressed  in  terms  of  b takes  the  more  familiar  form 

^ s - |G|  sin  6 
dt  ,, 

0 

(with  6m  0 , as  previously  explained,  and  gg  ° <*  lOj  sin  0 ) . 

As  is  veil  known,  (3*6)  and  (^.7)  <^on  be  integrated  by  quadratures 
in  terms  of  elemsntcxy  functions.  This  is  effected  by  eliminating  t 
and  using  p = UQg2  as  a nev  variable.  Thus  ve  get 

dp/dgg  = -f22PVs2^(/G/^  ” 

in  which  the  variables  are  separable.  Of  course,  the  actual  integration 
of  (^.6)  and  (5.7)  by  quadratures  in  the  Indicated  manner  aay  present 
practical  difficulties.  Approxliimte  solutions  ace  for  the  present  purposes 
quite  good  enough. 

Having  once  obtained  the  zerolib  approximation  for  u^  as  a functicm  of 

t,  the  zero^  approximation  for  COq  is  obtained  by  integrating  (5*2).  Wbjn 
this  has  been  accon5)li8hed,  so  that  both  u,.  and  a>  as  well  as  g,  or  rather 
their  zero^  approximations,  are  known  functions  of  t,  we  find  that  equations 
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(5.4)  iind  (5*5)  take  the  form  (4.1),  with  H(t)  = 0 and  L(t)  = g.  Good 
zerowJ  approximations  for  *vand  u can  therefore  be  obtained  from  formulas 

(4.4) .  In  making  this  statement  we  use  the  intuitively  obvious  fact 
that  both  Uq  and  have  con5)aratively  small  time  rates  of  change.  In 


other  words  our  zero"^  approximation  is  best  ^rtien  du^/dt  and  dcj  /dt 

0 0 . . • 

are  smallest . TSiis  will  insure  the  smallness  of  the  quantities  ^ 

of  the  preceding  section. 


6.  ^ the  most  favorable  metric  for  the  ballistic  equations.  The 
sensitivity  of  the  motion  of  the  shell  to  slight  variations  in  the  initial 
conditions  was  the  original  subject  of  this  research  project.  This  is 
closely  connected,  if  not  identical,  with  the  question  of  stability.  )ut 
so  far  the  studies  of  this  question  have  been  limited  to  systems  corrc- 
spending  to  equations  (5.4)  and  (5.5)  in  their  application  to  the  finding 
oi*  approximations.  In  such  applications,  the  Uq,  and  g are  regarded 

as  known.  In  fact  it  is  sometimes  assumed  in  effect  that  they  are  even 
constant.  The  condition  for  stability,  as  given  by  Nielsen  and  Synge5 
for  example,  is  then  to  the  effect  that  the  determinantal  equation 


(6.1) 


OL-T  P j 

y S-fj 


in  (T  should  be  such  that  the  real  part  of  each  root  should  be  negative 
(or  possibly  zero).  Here,  in  conformity  with  our  previous  notation,  we 
define  oc,  /?,  y , S'  as  follows : 


cx  = + i(ag^  h)6*.^ 

(6.2)  ^-^11^0  + 1^22’^0 

^ = *^11^0  ^(^22 
S z di2Uo  + i<i21^0‘ 

so  that  equations  (5.4)  and  (5.5)  appear  as  (4.1).  Such  a theory  of 
stability  does  indeed  take  into  account  the  part  of  the  mathematical 
theory  \Aich  one  intuitively  feels  is  the  most  crucial.  Nevertheless 
it  focuses  exclusive  attention  on  (5.4)  and  (5*5)  and  ignores  the  fact 
y ) S , and  g are  also  subject  to  disturbances,  since  the  two 

equations  in  question  are  really  just  a part  of  the  larger  system 
consisting  of  equations  (5.I)  - (5.5). 


■-.■V 
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In  order  to  handle  the  complete  system,  the  author  devised  a metric 
theory  (cf . MP-^  and  MPg)  which  may  he  applied  to  any  system  of  differential 

equations  of  the  form, 

dx^/dt  ^ tj  i ■ 1,  *•*,  n. 


In  this  theory  we  introduce  a Riemann  metric  (in  general  varying  with  the 
time  t)  for  the  n dimensional  space  of  the  x's.  An  inequality  is  then 
developed  for  the  "distance  between  two  solutions"  at  time  t in  terms  of 
the  distance  at  time  0.  From  this  inequality  rigorous  estimates  can  be 
read  off  to  answer  at  least  some  luestlons  about  the  way  in  which  the 
solution  depends  on  initial  conditions.  Results  are  certainly  obtained 
no  matter  how  the  metric  is  chosen  but,  unless  the  choice  is  a fortunate 
one,  the  results  may  turn  out  to  be  too  trivial  to  be  of  practical  value. 

A minimum  desideratum  would  be  to  choose  the  metric  in  such  a way  that  it 
would  be  sensitive  to  situations  causing  stability  or  instability  with 
regard  to  the  simpj,lfied  system  consisting  of  (5.4)  and  (5.5). 

For  bhis  reason  a preliminary  investigation  was  carried  out  to  choose 
a favorable  metric  for  the  study  of  (4.1).  It  my  be  noted  in  this  con- 
nection that  the  choice  is  not  influenced  by  the  presence  of  the  terms 
H(t)  and  L(t)  in  (4,1).  Biis  is  because  the  difference  between  two 
solutions  of  the  linear  non-homogeneous  equations  always  satisfies  the 
homogeneous  equations,  a theory  of  choosing  a metric  for  a linear  system 
is  indicated  in  Section  5 of  MPg.  When  this  theory  is  applied  to  (4.1),  we 

are  lad  to  the  following  quadratic  form  for  the  Riemaim  metric  of  the 
four  real  dimensional  space  of  iV=  i and  u ■ u^  + iUg. 

ds^^  - /dW“  + (2r  )"^(J -S)du/du  + (2  f )"^(  J - <x  )di?du  + 

(I  S'  -c^i^  + l(F-c<f  +4/?v|)  }2rl'^  /duf. 


where  iv,  u,  «,  y , etc.,  are  used  to  indicate  the  conjugate  Imaglnaries  of 

w,  u,  « , Y , etc.  In  another  paper^  this  metric  was  actually  used  to 
obtain  inequalities  for  the  system  (4.1)  and  these  inequalities  were 
shown  to  be  the  best  possible  ones  of  their  type.  Hence  we  are  confident 
th,.t  in  setting  up  a favorable  metric  for  the  entire  system  d h/  and  du 

should  enter  as  in  ds^^  or  (what  leads  to  the  same  result)  dsj|^2  multiplied 
by  on  arbitrary  positive  function. 

As  already  indicated,  the  rest  of  the  system  is  not  so  sensitive  to 
disturbing  phenomena  and  hence  there  seems  to  be  no  particular  reason 
for  using  anything  other  than  the  simplest  choice  of  a Euclidean  metric: 
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2 P 2 ? 2 

dSg  » du^  <+  dg  + dg  + 


It  may  be  noted  in  this  connection  that  the  non  linear  system  of 
Section  6 of  is  concerned  with  the  corresponding  problem  of  particle 

exterior  ballistics,  in  which  x and  y denote  two  coa5)onents  of  the  velocity 
of  the  particle.  Here  the  use  of  a Euclidean  metric  is  extremely  satisfactory. 

Ihus,  fer  the  co«q;>lete  system  one  would  exp*‘^t  satlsf action  from  a 
metric  of  the  form, 

ds2  e J2ii|  ^ds^2  ^ jds^. 

Here  dsj^^  is  multiplied  by  / 2Y|  ^ in  order  to  avoid  fractions  and  J is  a 

positive  constant  to  be  chosen  at  pleasure.  Considerable  effort  has  been 
spent  to  determine  an  efficient  value  for  J.  Since  these  efforts  have 
been  xiravalllns,  one  way  suppose  that  different  values  may  be  desired  for 
dlfferexit  aprjlicatlons . 

7*  Pie  fundamental  quadratic  tora»  We  turn  now  to  the  fundamental  result 
of  MP]_  lAlch  implies  that,  if^t)  demotes  the  distance  between  solutions 
at  time  t of  the  system. 


(7.1) 


dxj/dt  . X^(x^,  x^),  i « 1, 


where  "distance"  is  taken  in  the  Riemannian  sense  with  respect  to  a 
specified  metric. 


(7.2) 


ds^  * g .dx.dx  , 

1,J  iJ  i J 


' . ? 


(7.3) 


where  ^ is  an  upper  bound  for  a certain  quadratic  form 

■‘i.j  ^ -'y 

mdar  the  eenditlen  ^ = X.  acre  the  ceefflcieats  of  the 

"fundamental"  quadratic  form  Q^aJ  are  given  by  the  following  formula." 

(7.5)  ^ 


1:9; 


' lur, 
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In  this  section  we  wish  to  produce  a table  for  the  coefficients  of 

this  quadratic  form  for  the  ballistic  equations  (5.I)  - (5*5)  with  respect 

to  the  metric  introduced  in  the  preceding  section.  First,  however,  we 
indicate  the  change  of  notation  to  harmonize  MP,  with  the  notation  of 
the  exterior  ballistic  problem,  and  then  write  out  explicitly  the  differ- 
ential equations  and  the  coefficients  of  the  chosen  metric  in  the  notation 
of  MP^. 

In  our  problem  the  n of  MP^  is  equal  to  8,  and  we  choose  x^^  = 

Xg  -£*>2,  X3  - u^,  X|^  » Ug,  X5  . Uq,  xg  * - S2,  X0  -o/q.  Then  the 

differential  equations  in  the  new  notation  are  as  follows: 

^ = a^^x^  - (a^  + -f  ^11^3X3  - ■>>22X4X5 

dt 

^ « X = a^gXgX^  + (a^  + h)x^XQ  + h^X4Xg  + b2sX3X^ 

(!Biese  two  equations  were  obtained  by  equating  real  and  pure  imaginary 
part  of  (^•h)). 

^3  3 

X . ICj^xe  - (C22  + DXgXj  + di2  X3X5  . agiXi,Xg  + Xg 


dx, 

4 


a 


CiiXgXg  . l)x^x^  + d^^X^^Xj  * dgj^XjXg  + x^ 


(ihese  two  equations  were  obtained  by  e,,uating  real  and  pure  imaginary 
parts  x)f  (5.5)). 


dx^ 


dt 


= x5  = 


: xjXg  - x^x,^  - xpi  + f^^x^ 


dx, 

dt 

dx, 

at 


6 , - f Ini^  2 2 1 

^ : X2  (/g|  - Xg  -x^  )•■' 

I = xT  i . (|G)2  - Xg2  - X,y2)i 


*=8  . x8  . 

5J-  - X“  : 
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We  similarly  write  down  the  explicit  expressions  for  the 

Sll“  ^0  ■*"  '*0  ^ * ^^‘^22"*’^^  *5  ^8 

®12  ■ ° 

S13  * 2 f <=22^'^12"^12^  '*'  *5*8 

®21  ' ®22  ' ®n'  ®23  * "®1*^*  ®2l}-  ' ®13 

®31  ' 813^  ®32  ' “8l4 

833  = (V“12^S^+  '' 


where 


■ <V“12)"  - '«'22''=22+^>'  ^2  = '•Vn  - 


634  ■ 643^  ’ ej_4>  842  = 63^3/  84^  0,  g44  ■ g33 

®55  " «66  • 877  ‘ ®88  = 

gij  = 0 if  at  least  one  of  the  subscripts  is  > 4 and  if  i / J. 

On  account  of  relations,  dX^^Xg  = ^X^/dx^,  5X®/ ^Xg  =-^X^/^x^, 

?xV  = ^X^/7>x^f  3xV^Xg  = - "^X^/dx^,  etc.,  which  are  readily 

verified  from  the  formulas  for  the  X's,  we  find  with  the  help  of  (7 >5) 
that 

*11  “ ^2*  ^12  " “‘^1"  %3  " %k^  ‘^13  " 

^l4  ■ "*l23^  I3I  “ ^42"  ^32  ■ "^1"  ^^34  * "‘^3* 
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These  relations  merely  IJBply  that  the  part  of  the  quadratic  form  q[a1 

involving  only  a Hermitian  form  in>^j^  -f  and 

A 2 + and  this  Is  just  what  we  might  have  e^qpected  from  the 

pedigree  of  ^[aJ^  I*  is  not  necessary  to  ccmqpute  either  or  q^]^, 

since  it  is  only  the  syanetrical  part  of  the  matrix  (q^^)  which  is 

essential.  Hence,  for  i,j  = 4,  there  are  only  six  quantities  to  be 
computed,  namely 

(T.5.1)  q^i  - qg2  • ^^(Cg2  + - x^xi^  - (|g|^  - Xg^  - 

+ 2(Cg2  -f-  1)^(^2.2  * * ^“ll^^^l2  ^2  * 

(7.5.2)  qj^2  '124"  " ^12)  + (°22+^)(%  " ^21  " " ^*4  " ®l) 

+ |^(^22  ^^12  ■*■  (®22  ^^^*^21  " ^21  ” **" 
^^11^°22  ^^21^^2  ” ^12^^°22  ] *5^*8 


+ " ^°11^^^21  ■ ®21  " 

^ i. 

\diere  we  have  set  R]^  * (|g|^  - xg^  - 

(7.5.3)  q3j^-  q4g  = ■ ^12^  + ^°22^^^^*^21  " ^21  “ ^))(V3  " ^1^4  " Hi)^8 

^ [^^22  “12  2aj^2)(‘=U^<^2  ‘ ^2^  ^ ^«22  ^ ^^(*^21  ‘ ^21  " 

- 2(cg2  + l)(d^2  - aj^)(agj^  + h)  + c^(di2  " ^5% 

+ “ll^^l  " ^21  ■ h)(d2i  + agj^  + h)xg3  .f  °U®2^8" 
where  we  have  set  Rg  = + (2^j^^2  " ^3^K%^  + ^2%^)- 
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(7.5.4)  « -qg3  • 2(Cgg+  l)(d^  - ‘ 

+ 2(c22  + 1)  [(ai2  - ‘i*)(f22  ^12^  “ ^^22^^22  ■*  ^>]  ^5^ 

-2  j^(d^  - 

(7.5.5)  <143^  ■ -132  • 2(Cg2  .»■  1)(<1l2  - a^KVs  " V4  " ®1^5 

^ (Cgg  -f  l)(di2  " *12^ ^^22  ■*"  *12 

+ (=22  + ^V?- 

(T-5.fi)  <33  • «i^  ■ i(*2=‘3  • Vk  ■ ®j.)  |^®(*12  ‘ *12^%  ■*■  J^j 

1 2 ^ ^2  ^ 

^ * ^22*5  ^ “i2^*8  ^2^5 

’ *12^ ^^*12  ^22)(<^  ■ ®12^  ” ^22^°22  ’*’ 

Dm  etheir  <hi  ere  reXatlyely  sliqple  to  calculate  frm  (7*5)  end  vUl 
IM  left  to  the  reader.  The  conplexities  of  feraulaa  (7.5*l)  (7.5*6) 

pM»va  It  eee«  unllkelj  that  any  slnyle  general  concluaions  can  he  dravn. 

The  Method  glrea  prlnarlly  a procedure  for  the  nuaerlcal  study  of  indlrldual 

caeea. 


8.  A verd  of  vanlng.  After  writing  the  above  aatcrlal,  it  gradually 
cane  to  the  author's  attention  that  the  elimination  of  by  means  of 

g^2  ^ g^2  ^ g^2  ^ ^^2  ^ I Qj^  introduces  a difficulty  in  the  study  of  a 

trajectory  near  its  vertex,  that  is,  where  the  axis  of  the  shell  is  hearly 
hoarlsootal.  This  is  because  the  derivative  of  the  function  with  respect 

to  Xg  cr  x_  is  very  large  ^dien  g_  or  R^  is  close  to  0.  Hence  the  method, 
as  wQrlced  6ut  above,  is  pertinent  largely  to  antlalmraft  fire  or  other 
types  of  trajectories  where  the  shell  axis  is  never  horizontal  in  the 
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laportant  pert  of  the  trajectory.  If  this  is  not  the  case,  the  theory 
aay  he  andlfled  In  one  of  tvo  vays.  SLther  eliminate  Instead  of 

(this  has  the  disadvantage  of  making  our  equations  acre  unsymetrlcal) 
or  else  refrain  fren  eliminating  any  of  the  g's  and  use  a system  of 
order  nine  Instead  of  one  of  order  el^.  m either  event,  the  modi> 
flcatiens  necessary  for  this  purpose  are  only  sll|^,  but  it  is  not 
possible  to  vork  out  the  detidls  here. 
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